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Preface 

This is the fifth technical report covering the research work 
under the project entitled "A Test Program to Determine the Mechanical 
Behavior of Solid Propellants". The work reported here particularly 
refers to experimental studies on a rubber-like material that has 
been subjected to a uniaxial compression. The strains on the surface 
of a cylindrical specimen uncLer given loads were measured via the 
moir£ method of strain analysis. The boundary conditions on the 
ends of the cylinders were varied by varying the friction between 
the ends of the cylinders and the compression plates. 
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Abstract 

Experiments were conducted by compressing, uniaxially, a series 
of cylinders of length to diameter ratio of' two. The boundary 
conditions at the ends of the cylinders were varied with respect to 
the amount of friction present. The cylinders were cast of Solithane 
113, a rubber-like material. The compressive stresses applied varied 
from 0 to 19*1 lb/in 2 , for stresses which were based on the original 
cross-sectional area. Strains were measured on the surfaces of the 
cylinders, which corresponded to ten different values of applied 
stress. Strain was measured in the axial and circumferential 
directions by using the moire? fringe method. A short series of 
tension experiments were also conducted in order to verify the 
accuracy of the moir€ fringe method as employed for this presentation. 

A theoretical analysis was performed. The classical elastic 
solution to the problem of uniaxial compression of a cylinder has been 
done and is presented here in detail. A short review of the finite 
deformation of a cylinder of Neo-Hookean material under uniaxial com- 
pression which produces bulging is given. 

The theoretical strains as predicted by the classical theory are 
compared to the experimentally determined strains. A comparison of 
the strains obtained experimentally with different boundary conditions 
is presented. Also, a comparison of the stress- strain curves for 
Solithane in compression and tension is given. 
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Cl TMTER 1 

INTRODUCTION 

1.1 Statement of the Problem 

In recent years much attention has been drawn to new classes of 
materials. One important class , the rubber-like material, is already 
in wide use. As an example, solid fuel propellants are placed in this 
category. Certain rubber-like materials are also used as vibration 
dampers. In order to construct reliable engineering structures, it is 
important to know just how a class of materials will deform under a 
given system of loads and with given boundary conditions. 

The subject of this presentation is the study of a rubber-like 
material subjected to a uniaxial compressive load. Experimentally, a 
right, circular cylinder with length to diameter ratio of two is con- 
sidered. The boundary conditions are with respect to the amount of 
friction present at the plane ends of the cylinder where the load ap- 
plication occurs. The axial and circumferential strains on the sur- 
face of the cylinder at its axial center are measured as the load 
varies from zero to sixty pounds, in increments of six pounds. The 
end boundary conditions are also varied from one experiment to the 
next, i,e. , the amount of friction between the end plates and the 
specimen ends is varied by the application of a lubricant. 

A relatively new method of strain measurement, the moiri method, 
is employed, Therefore, tension experiments designed to measure modu- 
lus were also conducted and the results, using the moire method, are 
compared to the results of previously conducted tension experiments 
where clip gages had been used for strain measurement. This comparison 

6 




2 


is made in order to verify the accuracy of the moir£ fringe method of 
strain measurement, as used here. 

An attempt to solve the problem, theoretically, of the large 
deformation under uniaxial compression, of a right, circular cylinder 
of incompressible material is made. It is assumed that there i s no 
slipping of the ends of the specimen with respect to the compression 
plates during an application of load. Due to mathematical intracta- 
bilities in the solution of the problem, no final solution is given 
in this discussion. 

1 . 2 Importance of the Investigation 

A study of .the aforementioned problem is necessitated by several 
things. As has been mentioned, many so-called rubber- like materials 
are now in existence and, conceivably, many more will be developed in 
the future. In order to design reliable structures with a minimum of 
material and cost, and having a maximum of safety, the material proper- 
ties of the design material must be known. 

With regard to this specific problem, compression is surely one 
of the important loading configurations to be considered in most de- 
sign problems. For instance, the solid fuel propellant in a rocket 
sometimes has the shape of a thick-walled cylindrical tube. The burn- 
ing of the fuel proceeds from the inner wall to the outer wall. With 
this burning, there exists high temperatures and pressures. These 
high temperatures and pressures exert compressive forces on the wall 
of the cylindrical tube. As can be seen by this example, compression 
is an important loading configuration and this importance warrants 
further study of the problem. 
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1 .3 Scope and Limitations 

The work presented here delves into one aspect of the problem of 
compression of materials. The geometrical shape of the body being 
considered, the manner in which the load is being applied, the 
boundary conditions being imposed, and the type of material being 
considered all contribute to the limitations of the scope of the 
problem, both experimentally and theoretically. 

A right, circular cylinder, with the length being twice the 
diameter, will be considered, both experimentally and theoretically. 
This length to diameter ratio is chosen in order to facilitate good 
bulging of the curved surface of the cylinder. A cylindrical shape, 
instead of a parallelepiped, is chosen solely for reasons of theo- 
retical simplicity. A parallelepiped would yield infinite stresses 
at its edges. 

A uniaxial compressive load is applied to the terminal ends of 
the cylinder in both the theoretical and experimental investigations. 
In the theoretical case, the total load is assumed t o be known, how- 
ever, the manner in which this load is distributed over the end sur- 
faces is unknown, except that at the edge of each surface the load 
must be zero. Also, the total load applied is considered in the ex- 
perimental investigation. The total load is varied from zero to sixty 
pounds and the strains are determined experimentally and theoretically 
at increments of six pounds. 

For the theoretical case, the boundary conditions assumed are as 
follows. A given total stress acts on the ends. Nd axial or cir- 
cumferential stresses exist on the curved surface of the cylinder. A 
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shear stress of unknown distribution acts on each end in order to con- 
strict the ends from axial movement as the load is applied. Also, the 
ends are constrained to remain plane when the load is applied. 

For the experimental case, one of the boundary conditions is 
varied fromthat in the theoretical case. A given total stress is 
applied to the ends as previously mentioned. Also, no axial or cir- 
cumferential stresses exist on the curved surface. The ends are con- 
strained to remain plane as the load is applied. But, the shear 
stress on the end is not always the same as in the theoretical case. 
Two cases are considered. The ends are constricted from axial move- 
ment in one case. In the other case, there is some amount of axial 
movement of the ends for a given load. 

In the theoretical case, the material is considered t o be Neo- 
Hookean, A Neo-Hookean material is one that obeys a type of stress- 
strain relationship that is shown to be a natural extension of the 
Hooke's law used in the study of small elastic deformations. The 
material is also assumed to be homogeneous, isotropic, and incom- 
pressible. 

Experimentally, the material considered is similar, in some 
aspects, to that assumed for the theoretical treatment, ma terial 

used.; Solithane II3» is homogeneous and isotropic. As to whether 
Solithane 113 is incompressible and a Neo-Hookean material is still 
unknown. 

Chapter Two deals with the theoretical investigation. Part One 
is a study of the classical elasticity solution to the problem of a 
cylinder of moderate length being compressed uniaxially with given 
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boundary conditions as derived by Filon (3). Section Two of Chapter 
Two is a short review of the finite deformation theory as presented by 
Rivlin (6). The basic equations of this theory are set forth and in- 
clude the stress versus deformation equations, the incompressibility 
condition, the equations of equilibrium, et al. Part Three of this 
chapter is an application of the finite deformation theory to the 
specific problem dealt with in this presentation. Part Four sets 
forth the results of the calculations of the surface strains as pre- 
dicted by the classical elasticity theory. 

The experimental investigation is described in Chapter Three. 

A chemical description of the material used in the experimental 
program is given in Section One of Chapter Three. Part Two of this 
chapter deals with the exact manner in which the specimens that were 
used in the experimental investigations were made. The procedure 
used in testing the specimens is described in Part Three of this 
chapter. Also in this part, one will find a description of the 
various apparatus used in the experimental program. Part Four de- 
scribes the theory behind the moir£ method of strain analysis. This 
description includes a short review of three of the most common 
methods of analyzing moir€ fringes along with a comparative analysis 
of the three different methods. 

In Chapter Four the results of the experimental investigations 
are set forth. Data for both tension and compression experiments are 
presented. A comparison between the experimental results in com- 
pression, with the proper boundary conditions, and the theoretical 
results as given by the classical theory of elasticity is made. All 
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prevalent results of the experimental data are discussed in detail, 
and conclusions are drawn. 
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CHAPTER 2 

THEORETICAL INVESTIGATION 


2.1 Solution of the Problem by Classical Elasticity Theory 

The following was taken from parts of a paper by L.W-.G. Filon (3). 
The reasons for reproducing it here are two fold. First, as presented 
here, the solution of the problem is in much greater detail than in the 
original paper. Also, the problem, as solved by Filon, was studied in 
order to obtain some insight into the problem of the compression of a 
cylinder in general and the boundary conditions of that problem in 
particular. 

The problem discussed is that of a cylinder of moderate length, 
when compared to the diameter, which is compressed between two rigid 
planes in such away that the terminal cross-sections are constrained 
to remain plane. In this section of the theoretical development, the 
general equations for the elastic case of a right, circular cylinder 
of finite length are obtained. The method adopted has been to obtain 
symmetrical solutions of the equations of elasticity in cylindrical 
coordinates and to express the typical term in the form (cos kz)xf(r) 
and (sin kz )xf (r ) . 

Let us first consider the differential equations of equilibrium 
in cylindrical coordinates as shown by Timoshenko (9). 
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3a, 3a a 

®z zz rz 
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9ct a , 9a, , 9a, 2 a , 

r<? + 1 <j)j> + <(> z + r<j> _ 

9r r 9<j> 9z r 


( 2 . 1 . 1 ) 


Writing the stresses in terms of displacements and employing Lamd's 
constants 


rr 


(A+ 2 „) M + X " + X 
v p dr r dz 


, du' .. dw . u 

a = A — ■{ A+2y ) z — — 

zz dr dz r 


a, , = 




U+ 21.1) f + A g + X * 


dz 


rz 


= P 


dw du 
dr dz 


<pz 


dv 

dz 


r<j> 


dv v 
dr r 


(2.1.2) 


Placing the stresses as written in equations (2.1.2) into the equili- 
brium equation s (2.1.1) , we obtain after rearranging and assuming that 
u, v, and w are independent of <j>, 


( A+2y ) ^ + ( X+2y ) d ~ © + „ ^ + ( X+p ) - 0 (.) 

dr 2 dz 2 drdz 


d fl d(rv)l 


dr I r dr 


d z v 

dz 2 


U+p) 


d 2 u ( 1 du) 


d 2 w 1_ dw) 

drdz r dzj 

+ P 

,2 r 
k dr z J 


+ (X+2y) 


dfw 

dz 2 


= 0 . 


(b) (2.1.3) 

(c) 


Let us define the operators 


and 


£ 2 E 


did 
dr r dr 


D 


d_ 

dz 


(2.1.4) 
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Differentiating (2.1. 3. a) witl} respect to z and (2.1. 3. c) with 
respect tor. 


(X+2y) 


d 

d 2 u 

i ( \ i Q , - \ _ 

— f-1 " 

+ „ * 

d 2 u' 


d 2 w ' 

dz 

dr 2 . 

+ (X+2y ) dz 

dr '•r' 

+ y & 

dz 2 

+ [x+]1 > dT 

drdz 


= 0 


and 

(x+y) 


d 

d 2 u ( 1 du 

+ „ £- 

d 2 w 1 dw 

+ ( \ 4- On ^ - 

d^ 

dr 

drdz r dz. 

y dr 

dr 2 r dr . 

+ (X+ 2y) ^ 

dz 2 


= 0 


+ (X+y) D 2 ~ = 0 
dz dr 


(X + y) ♦ 


dw 

dr 


= 0 


Rearranging and substituting equations (2.1.4). 

[( A+ 2y ) £ 2 + y D 2 

y £ 2 + (A+ 21.1) D 2 
and from equation (2.1.3.b) 

U 2 + D 2 ) v = 0. 

du 

From equations (2.1.5) and (2.1.6) eliminate and rearrange 


(2.1.5) 

( 2 . 1 . 6 ) 


(2.1.7) 


-(X+y ) 2 £ 2 D 2 


dw 

dr 


(X+2y)y a h + (X+2y)yD 4 + y 2 £ 2 D 2 + (X+2y ) 2 £ 2 D 2 


dw 

dr 


[[(A+ 


2y)y]£ 4 + 2y(X+ 2y)£ 2 D 2 + (X+ 2y)yD 4 


dw 

dr 


0 


[£ 4 + 2£ 2 D 2 + D 4 ] ^ = 0 


( £ 2 + D 2 ) 2 ^ = 0 = (£ 2 + D 2 ) 2 y l 
Similarly, by eliminating ~ 

(£ 2 + D 2 ) 2 jg = 0 = (£ 2 + D 2 ) 2 
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The solution of this differential equation determines the elastic 
equilibrium of the circular cylinder under any symmetrical system of 
stresses. 

The above differential equation can be solved by the separation 
of variables method. Consider first the equation (f 2 + D 2 )y = 0. Let 


y 3 = R 1 (r)Z 1 (z) 


U 2 + D 2 )R 1 (r)Z 1 (z) = 0 


1 d 


a 2 


2+ r r 2 da 2 


dr 
d 2 R x 


R 1 (r)Z 1 (z) = 0 
d 2 Z, 


.. dRj R^Z^ ll zj ^ 

z £ — _ + R 

l r dr t r, 

dz 


Dividing by Z to separate the variables. 


d 2 R x 


dRi ! ! d 2 Zi 


1 dr 2 R l r dr r 2 Z 1 dz 2 


= 0 


Rearranging 


1 + i i 

R, i , 2 R x r dr ” r2 


d 2 Z, 

1 = + it 2 


T dr" 


Z 1 dz 2 


Separating the equation 

i 

E i dr 2 v 


dR 

dr 


1 1 


= + k 2 


1 




dz 2 


- - k 2 . 
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Rearranging 


d 2 R 1 dR 1 

, o + r dr 
dr z 


i— + k2 


R = 0 
1 


( 2 . 1 . 8 ) 


whose solution is 

R 1 = AI^Gtr) + BK^kr) 


d 2 Z 


dz 


- ■+ k2 z - o 
2 1 


( 2 . 1 . 9 ) 


Since y 3 is the solution of the equation 


U 2 + D 2 ) y = 0 , 


pick y so that it is a solution of the equation 


(£ 2 + D 2 )y = y 3 . 


Multiplying both sides by (£ 2 + D 2 ) 


( £ 2 + D 2 ) y = y 3 U 2 + D 2 ) = 0 


So that, y is also a solution of the equation 
4 


U 2 + D 2 ) y = 0 . 


Solve the equation 


U 2 + D 2 ) y 4 = y 3 


Let 


y 3 = R 3 (r)Z 3 (z) and y 4 = R 4 (r)Z 4 (z). 


(£ 2 + D 2 )R Z = R Z 
4 4 3 3 
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Rearrange 


£ 2 R 


R. 


» ¥j 

\ ' E A 


( 2 . 1 . 10 ) 


This is a function of r plus a function of z equal to a product 
function of r and z. If “equation (2. 1 . 10) is to be satisfied identi- 
cally, the product function must be a function of either r or z only. 
(Case i ) 

Let = aZg, where a = constant. 

R 3 Z 3 

22 + 1,2 = v ^ 

L <* 2 + d2 > -is. 


From equation (2.1.9) 


and. 


(p 2 + k 2 ) Z 1 = 0 
D 2 = - k' 2 . 


Therefore, 


-K 2 ) =fE 3 


Multiply both sides by (£ 2 ~ k 2 ). 


R {l 2 ~ k 2 ) 2 = £ R (SL 2 - k 2 ) = 0 
4 a 3 


Therefore, R, is a solution of 


U 2 " k 2 ) R. = 0 


From equation (2.1.4) 


£ 2 = 


d 2 1 d_ 1_ 

dr 2 r ^ “ r 2 
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a 2 

dr 2 


+ 


i a_ 
r dr 




The method of solution of this equation is .presented in Appendix A. 

The solutions of this equation are 

I^kr), K^kr), I^kr), —■ K^kr). 

Since 

££ I]_(kr) = rl ^ (kr) = rI Q (kr) - t H(kr) 

and 

J ^(kr) = rK Q (kr) - £ K^kr), 

the four independent integrals are 

I^kr), KjOcr), rI Q (kr), rK Q (kr). 

(Case ii ) 

Let R = bR 0 . 

H 0 

From equation (2.1.9) s ' 

(D 2 + k 2 ) Z l = 0 (2.1.11) 

d 2 Z 

I + k 2 Z = 0 

dz 2 

D 2 = - k 2 
D = ± ik. 

Therefore, 

Zj = B^sinkz + B^coskz (2. 1.11. a) 

where B^ and B 2 are constants. 
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Let = bR 3 in equation (2.1.10) where b equals a constant. 


U 2 + d 2 )z 4 = z 3 /b 


From equation (2.1.8) 


U 2 " k 2 )R L = 0 or l 2 = k 2 . 


Or, substituting in the above equation 


(D 2 + k' 2 )Z 4 = Z 3 /b. 


Multiply both sides by (D 2 + k 2 ). 

2 

(D 2 + k 2 ) Z, = £ (D 2 + k 2 )Z $ = 0 

Assume the solution 


Z 4 = A coskz. 


Substitute in the following equation where ~ = Bcoskz from the above 
solution (2.1.11. a;. 

(D 2 + K 2 )Z = Bcoskz 
1 4 


Since 


and 


Z, = Acoskz, 

4 ’ 

D 2 Z, = “ k 2 Acoskz, 
4 


Z 4 (- K 2 Ac.oskz + k 2 Acoskz) = Bcoskz 


cannot be solved this way, therefore use the fact that 

ikz . -ikz 
, e e 

coskz = n 


a 



in the equation 


Z„ r 


D 2 + k 2 


coskz. 


Z„ I 


1 

ikz , -ikz 
e + e 

D 2 + k 2 

2 


i ikz -ikz 

Z = x ^ x 

4 D 2 + ^2 ^ D 2 + k2 2 


ikz 


Z, = 


1 x ® 


4 ikz 


4 (D+ik) (b-ik) 2 (D+ik)~(D-ik) 


Z, = 


Z„ = 


1 

ikz 

e 

1 

-ikz 

e 

D-ik " 4ki 

~ D+ik 

* 4ki 

k£r 

ikz 

eakz 

-ikz 

e 


D-ik “ 

D+ik 



Z 4 4ki 


^ikz -ikz 


Z 4 = ^ zsink2 


Similarly, it can be shown that 


Z, = — z coskz. 
4 2k 


Therefore, the possible sets of product functions that satisfy 
2 
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y = A cos (kz + a) I^kr) 

y = B cos (kz + g/ K^kr) 

y = C cos (kz + y) r I (kr) 

( 2 . 1 . 12 ) 

y = D cos (kz + 6) r K Q (kr) 

y = Ezcos (kz + E) I (kr) 

y = Fzcos (kz +9) K^(kr) 

where a, 3 S Yj 6, e, 0 are constants. 

Let us consider the general case of the following system of stresses: 

(i ) a rr /u = a §^ ven even fctn of z over the curved surface r = a. 

(ii_) CT r2 /u = a given °dd fctn of z over the curved surface r = a. 

(ill.) cr #0 over the plane ends z = * c. 
z z 

2 

Since du/dz and dw/dr are both solutions of (£ 2 + D 2 ) y = 0 we may 
express them by a series of terms as follows: 


dn- =1 { AjCos (kz + o^Jl^kr) + C x cos(kz + Y 1 )rI Q (kr) (2.1.13) 
+ EjZ cos(kz + Sjjl^kr)! 


err = l ^ A 2 cos(kz + a^I^kr) + C 2 cos(kz + y2)rI Q (kr) (2.1.14) 
+ E 2 z cos(kz + e 2 )l 1 (kr) j . 

No K-functions have been ihtroduced, as they lead to infinite terms at 
r = 0. 

It can be reasoned that, for a cylinder u should be an even function 
of z and that w should be an odd function of z (u is symmetric about 
the z-axis and w is skew- symmetric about the z-axis). 
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Let 

<*1 = a 2 = Y t = y 2 = -tt/ 2 ; e x * e 2 ** 0 

The reason for this will be seen upon integration of equations (2.1,13) 
and (2,1,14) as follows, 

= l ^ A 1 cos(kz - ir/S)l 1 (kr) + CjCosUz - n/2)rl 0 (kr) 

+ E^z cos(kz.)lj (kr) j 

dw = l { A 2 cos(kz - 11 / 2 ) 1 ^ kr) + C 2 cos(kz “ ir/2)rl 0 (kr) 
dr 

+ E 2 z cos(kz)lj(kr) ^ 

Expanding the trigonometric functions in squat io;njs (2.1.13) and 

(2.1.14) 

^ ^ Ajfcoskz cas(- ir/2) + sinkz sin(Tr/2) ]I ^ (kr) + 

C [coskz oos(- tr/2 ) + sinkz sin(ir/2) ]rl^(kr) + 

E^z cos(kz)l ;L (kr) 


- £ ^A 2 [coskz cos(—nr/2) + sinkz sin(ir/2) ]I 1 (kr) + 

C 2 [coskz cos(- tt/2) + sinkz sin(Tr/2)]rI Q (kr) + 

E 2 z cos(kz)l 1 (kr) I 

.Reducing 

On = M + A 1 sin(kz)l 1 (kr) + C 1 sin(kz)rI Q (kr) + E 1 zcos(kz)l 1 (kr) 
dz 1 


(2,1,15) 
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= A 2 sin(kz )l x (kr) + C 2 sin(kz)rl 0 (kr) + E 2 zcos(kz)l 1 (kr)^ 


(2.1.16) 


Integrating (2.1.15) with respect to z and (2.1.16) with respect to r 


u 




— cos(kz)l 1 (kr) - ^ cos(kz)rl 0 (kr) + 


E 


z sin(kz) + 


cos (kz )] 


I x (kr) j + x(r) (2.1.17) 


w 


sin(kz)l Q (kr) + ^2" sin(kz)rl 1 (kr) 

E 'i 

+ z cos(kz)l Q (kr) >■ + 0(z) 


(2.1.18) 


To evaluate the constants, substitute the proper derivatives of 

( 2 . 1 . 17 ) and ( 2 . 1 . 18 ) into equations (2.1.1) and (2,1,3). 

(X+2y) Jj j -Ajk cos(kz)l 1 "(kr) ~ C^k cos(kz)rl 0 "(kr) “ 

C cos(kz)l ' (kr) - CjCostkz)! '(kr) + 


E z sin(kz) + - °- s ^ — — kl "(kr) + A cos(kz) - — I,(kr) 

1 i it 


1 


- A 1 cos(kz) r 1 I 1 ' (kr) ~ C 1 cos(kz)lQ l (kr) 


+ E.. 


z sin(kz) + 


1 r~ 2 

r i I 1 ' (kr )- ~~ I x (kr) 


+ d 2 X.(r) 

dr 2 


+ x(p) 
r dr, 0 
r z 


+ y 


+ U+y) 


£ ^ + kA^os (kz (kr ) 

+ kC 1 cos(kz)rI 0 (|cr) + E^-k z sin(kz) + cos(kz) I 1 (kr)| 
J ^ A 2 cos (kz )kI Q 1 (kr) + C 2 cos(kz)(kr)l i ’ (kr) + 


i 
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+ C 2 cos(kz )l 1 (kr ) " E 2 z sin(kz )kI Q 1 (kr ) + 


E^cos 


(kz)l 0 ‘ (kr)j 


(2. 1.3. aj 


(A+p) A t sin(kz )kl 1 1 (kr) + C 1 sin(kz)(rk)l 0 ' (kr) 


+ C,sin(kz)l n (kr) + E^I^kr) 


1 A 

(kz)cos(kz) + 1 sin(kz )l . (kr ) 

j r i 


+ C 1 sin(kz)l 0 (kr) + ^ (z cos(kz))l 1 (kr)| J + p 
£ | kA 2 sin(kz)l 0 " (kr) + C 2 ksin(kz)rl l "(kr) + CgSir^kz)]^ * (kr) 


A 


+ C 2 sin(kz)l 1 1 (kr) + E 2 z cos(kz)kl 0 "(kr) + sin(kz)l 0 ' (kr) 

G E ' 'll 

+ C 2 sin(kz)l 1 ' (kr) + sin(kz)l 1 (kr ) + ^ z cos(kz)l 0 ' (kr )>l 


+ ( A+2p) ^ -A 2 sin(kz)kl 0 (kr) - CgrksirKkz)]^ (kr) 

- E 2 z cos(kz)kl 0 (kr) - E 2 sin(kz)l 0 (kr) - E 2 sin(kz )l Q (kr ) 


d 2 9 (z) 


dz' 


}] =0 - 


(2.1.3.c) 


Combine like terms in ( 2 . 1 . 3 , a ) and ( 2 . 1 . 3 . c ) 


' r ' 2 


-1. 


( A+2p ) I I|a ]l cos (kz ) I + - — I : ( kr ) - kl^kr) - r" I 1 '(kr) 

- C 1 cos(kz)|(kr)l ( , "(kr) + 3I 1 (kr) + Ej z sin(kz) 
cos (kz ) 


+ P 
+ E 


f " 2 -ill 

k I t " (kr)- ^ — 1^ (kr) + r ^'(kr) j + i 2 x(r) 

1 | +k A 1 cos(kz)l 1 (kr) + (rk) C 1 cos(kz)l 0 (kr) 

-(kz)sin(kz) + cos(kz) 


I r (kr)H + ( A+p ) x 


i 
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£ | A^cos (kz Jklj. (kr) + C 2 cos(kz) (krjl^kr) 


+ I x (kr) 
■ 0 


- E, 


^ ,V 

z sin(kz)kt 1 (kr) - cosCkzJl^krjli 




^j+A^intkz) kl 1 '(kr)+r 1 I 1 (kr) +C,sin(kz) 
l. n r * 


U+y) 


(krjl^kr) + 2I 0 (kr)j + E^oslkz) l(kz)l ■ 1 (kr) 


+ (r 1 z)l 1 (kr) 




r r 

)}] 

+ y 

J ( A 2 sin(kz) . 


- 1 . 


- 1 . 


+ C 2 sin(kz) (kr)X 1 "(kr) + 31^ (kr) + (kr) Ijfkr) 


E 2 cos (kz ) 


(kz)l 0 "(kr) + (r - z) I 0 '(kr) 


(X+2y) 


1 


{" A 2 sin(kz)kl 0 (kr) - C 2 sin(kz) (kr )l i (kr ) - E 2 
(( kz )cas (kz)l 0 ( kfj + 2sin(kz)l 0 (kr)\ | + D 2 e(z) = 0 


In (2. 1.3. a) and (2.1, 3. c) the recurrence relation I ^ (kr) = I 
has been used. 

Jt can be shown that we can reduce (2.1,3, a) and (2,1, 3. c 
the following four equations. 


£ 2 X (r). - 0 
D%(z) = 0 

(A 1 -A 2 )(X+y)k 2 + 2k{G 1 (X+2y)-y ^ - (X+y)E 2 } = 0 
(A 1 -A 2 )(X+y)k 2 + 2k{G 1 (X+y) +• yC 2 - U+2y)E 2 } = 0 

If we let C, = C 2 = C and E, = E 2 = E, then equations (2.1.21) 
(2.1.22) become one equation, as shown 


i 


t (kr) 

) to 

(2.1.19) 

( 2 . 1 . 20 ) 
( 2 . 1 . 21 ) 
( 2 . 1 . 22 ) 

and 
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(A x - A 2 )(X+y)J£ 2 + 2k{c(A+2y) " (A+2y)E} = 0 


(A x - A 2 )(A+y)k + 2 (l+2y)(C-E) = 0 . 

By integrating equations (2.1,19) and (2.1.20) we obtain 

x(r) = uor and 0(z) = woz . 

From equations (2.1.2) 


, , , „ v du . , u . du. 
0 = (A+2y) j* - + A — +A , 

rr v dr r dz 


a rr = ( A+2y ) l | ~A 1 cos(kz)l 1 * (kr) - C cos(kz)(r)l 0 ' (kr) 

ill 

- C cos(kz) (k -1 )l (kr) + E z sin(kz) +k cos(kz) I ' (kr) f 

0 L J 1 J 

+ dx(r)/dr + X M -A 1 .cos(kz ) (kr) -1 ^ (kr) ” C cos(kz)k _1 I0 (kr) 

— 1 ^—1 1-.1 
+ E Z sin(kz) + k cos(kz)^(kr) I^(kr) J + r x( r ) 

+ X ^|A 2 cos(kz)l 0 (krl + C cos(kz) (r)l 1 (kr) - E z sin(kz) 

I Q (kr) “ k -1 cos(kz)l 0 (kr) j + d0(z)/dz 

ct = (2A+2y)u + A w + 7 -2y C cos(kz)(r)l (kr) -2(X+y) 

rr ° o L 1 

C cos(kz)(k _1 )l Q (kr) + A Ez sin(kz)l 1 ' (kr) + AE k^costkz)^ ' (kr) 

+ 2LZEz sin(kz )I^ ' (kr) + 2y E k~’' 1 "cos(kz)l 1 ’ (kr) + 

+ A Ez sin(kz)(kr) -1 I 1 (kr) 

+ A E k~ 1 cos(kz)(kr)~ 1 I 1 (kr) ~ A Ez sin(kz)l 0 (kr) + 

+ A E k -1 cos(kz)l Q (kr) 

r \ 

- 2X A 1 cos(kz)l Q (kr) + 2y -AjCos (kz )l x ' (kr) + X A 2 cos(kz)l Q (kr) 

\ J 


j 
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This reduces to 


or r = 2(X+y)u Q + A w Q +1 - \ [ (2A+3u )A X + pA 2 | I Q (kr)coskz 

I,0tr)O 


+ 2pi A x - E/k 


'i Ij(kr) - 1 . 

— ; coskz + Ez sinkz I (kr) - — - 1 > 

kr ( 0 kr. JJ 


- Crl (kr )coskz 


Likewise, 


a = 2Au + (A+2u)w + 


zz 




(A+2u)A 2 “ A A, - A 


2c 


r 

2(A+|i) I (kr)coskz + 2\i ■ Crl. (kr )coskz 
u *• 


" E I (kr)z sinkzl 


and 


rz 


= y £ + A 2 )I 1 (kr)sinkz + 2CrI 0 (kr) sinkz 

+ 2EI 1 (kr)z coskz | 


(2.1.24) 


(2.1.25) 


(2.1.26) 


For the displacements let E 1 = E 2 = 0 and C, = C, = C and substitute 
x(r) = u Q r and 8 ( z ) = w q z into equations (2. l.lj) and (2.1.18). Then 
with kr = p 

u = ur + ^ ^ I x (p) - £ rl 0 (p) J coskz 

( 2.1 IT, a) 

W = w Q Z + l | ^ I 0 (p) + £ rX x (p ) 1 sinkz. 

, (2.1.18 ,a) 


Let us consider our cylinder so that the coordinate system may be 
defined and that the boundary conditions may be established. The co- 
ordinate system is as defined in Figure 1. 

3 
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The boundary conditions are as follows: 


(i) A total pressure ira 2 Q exists over the plane ends, the 

distribution of this pressure being, however, unknown. 

(ii) w = constant whep z = ± c. 

(iii) u = 0 when r = a, z = * c. 

(iv) o = 0 when r = a, 

a =0 tffeep r = a. 
rz 


Filon has shown that a solution of the form shown in equation 
(2.1,17. a) and (2.1.18. a) plus a finite power series in r and z will 
satisfy the boundary conditions, incompressibility condition, and 
equilibrium equations which have been specified. This assumed 
solution is given below. 


r A 

u = l {" iT V kr) 


£ ^ u^r^ Drz^ 

k rl n (kr) j coskz + Uqt + 3 + ~ 


Er 3 z 


3,2 


Frz L 


( 2 . 1 . 27 ) 


and 


n W Z 3 W Z 5 

= H Iq ( kr) + k- rl 1 (kr) V sinkz + w q z H 3 — + 


w 


O + 

D 1 r 2 z E 1 r 2 z 3 . F 1 r 4 z 


( 2 . 1 . 28 ) 


These sums of series have been restricted by the conditions that u 
must be odd in r and even in z, and that w must be even in r and odd 
in z . 

The method of solving for the unknown constants, that was em- 
ployed by Filon, cannot be used in the problem that is ultimately 
considered here. Filon was dealing with linear equations and the 


1 
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equations to be considered here are non-linear because this study in- 
volves a rubber-like material which has the characteristic of large 
deformations. For this reason Filon's method of solution will not be 
given here because it would not contribute anything toward a better 
understanding of the problem, Instead, the constants of equations 
( 2 . 1 , 27 ) and ( 2 . 1 , 28 ), as found by Filon, will be presented. 


u 

o 


Qg(2.Y-l) 

y [h + (2y-l)f + 2g(4y-l 


IT 


3 

u i = 8 


(l-y )a 2 - — c 2 


E 


u o = - 2 I (i-y)e 


D = (l-y)c 2 " a 2 | E 

F -Q(2_Y-1_) 

' p[h + (2y-l)f + 2g(4y-l)] 


F = - | (l+y)E 

Qh(l-y) + 2gy 

° y[h + (2y-l)f + 2g(4y-l)] 



(2.1.29) 


W 2 = 
D 1 


| ye 


4yc* 


E 


E 


1 _ 


yE 


F 1 = 0 
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where 


k 



n=l } 2 , 3 j » • • * 


(2.1.30) 


f = (2 Y -1) f- 


+ ~ a^c 


2^2. 


lj (2y + l)c 4 


g 


= (l- Y -^)a 4 - \ a 2 c 2 + -j 


1 f 


1+ Y - 3J [2y+1] 


h = 


6- l 


l4 a a 4 . ,\ . 2a 2 c 2 , Y 
— c 4 Y - 3J (2y+l) + ~ 3 — 

2l+ [ ( 2 y+ 1 ) a + £ (8 y+1)]I^ - (8y+4)I q I 1 


1 a 


Y<x 2 I 2 - I 2 (l+Y« 2 ) 
0 1 


(2.1.31) 


A+U 

X+2\i 


a = ka 


mra 

c 


In the above. 


= I 


2l+ [(2Y+l)a + — (8 y+1)]I 2 - (Sy+4)!^ 


1 a 


Y« 2 ! 2 - I 2 (l+Y<* 2 ) 


(2.1.32) 


We also have , 




. . <-p.) n+1 

Y E4a 2 

[(2 Y +l)a + 4 -^- + --]l^ - 4 I o I i (2 Y +1) 

, (2.1.33) 

C n = 2 


Y<x 2 I 2 - (1+yo 2 )I 2 

A 

2n - 2 

i6ac 3 
n 3 ir 3 

2aJ ' 

( 2t+ 1) B + C n ^ ^ 

(2.1.34) 


and 



1.53 



(2.1.35) 
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2.2 Review of Finite Deformation Theory 

The basic equations of the finite deformation theory were derived 
by E. and F. Cosserat in 1896. The following equations were set forth 
by Rivlin (6) in 1948. They are set forth here as a foundation upon 
which the theoretical portion of this present work rests. Rivlin has 
introduced the concept of a neo-Hookean material. This material is 
assumed t o be incompressible and isotropic in its undeformed state 
and is capable of large, elastic deformations. Rivlin has defined the 
stress- strain relationships for this neo-Hookean material and using 
these relationships has developed a mathematical theory for trite; ma- 
terial. The principal results obtained have been expressed with 
reference to a cylindrical polar coordinate system (r,8,z) and have 
been set forth in the following. 

The strain components are defined as being at the point which is 
at (r,8,z) in the undeformed state, as shown in Figure 1. 
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The incompressibility condition, using a cylindrical polar co- 
ordinate system, is 


-tls- 

1 

r 

( 

7 

7. 

1- T ■ 

8u. 

9z 

8v 


_i_ 

U 1 8v 

t r 30 

8v 

8r 

-L 

i 

8z 



1 

8w 

! + ^ 

8r 


r 

80 

8z 


( 2 . 2 . 2 ) 


Since, in the deformation, the point (r,0,z) moves to (r+s, 0+<j>, 


z+w) , 

u = (r+s) cos$—r, and v = (r+s)sin<j>. 


(2.2.3) 


Therefore, it has been shown that the stresses in terms of displace- 
ments can be written as. 
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The equations of motion are 
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where 
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(2.2.6) 


In the equations of motion, p Q is the density of the material and 
Rj, 0, and are the componehts of any body forces on the material in 
the r s 8, and z directions respectively. 

If the surface of the body considered has, in the undeformed 
state, the shape of a cylinder with the z-axis as the axis of symmetry. 
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the boundary conditions will be different on the plane and curved sur- 
faces. Rivlin has shown the boundary conditions to be, for the curved 
and plane surfaces respectively, 
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These equations can be simplified greatly, viz. If, as in the 
problem considered here, the deformation possesses cylindrical 
symmetry about the z-axis!, then u>VjW, and p are independent of 0, 
i . e . , 

8u 3v _ 3w 3p _ Q 
30-38 30 = a0 ~ U * 

Introducing this into equation (2.2.1) for the strains and with the 
equation v = 0 for all points (r,S ? z) because no rotation about the 
z-axis occurs, 
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The incompressibility condition becomes , 
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The equations of equilibrium become, with R, = 0 = = 0 and 
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The boundary conditions become, for the surface parallel to the z-axis 
and cylindrical in shape, 
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(2.2.14) 


and for the plane surfaces which are normal to the z-axis, 
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(2.2.16) 


For the curved surface and Z have been taken as ze;ro because for 
this specific problem that is the case. 


2.3 Finite Deformation Theory Applied to the Problem 

In the previous section the equations describing the stresses, 
strains, etc. of a right circular cylinder under an arbitrary system 
of loads have been presented. In this section an attempt to solve 
these equations for prescribed boundary conditions and a prescribed 
system of loading has been attempted. The body considered, along 
with the coordinate system used in the solution, is shown in Figure 
1. In this derivation 2c = L. 

The displacement boundary conditions are 
(i. ) w(r. If) = K = constant, 

C i i - ) u(o,z) = 0, 
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and 

(iii. ) u(a,|9 = 0, 

The third condition states . that there i s no displacement in the radial 
direction, at the edges (r=a) of the ends of the cylinder. 

The stress boundary conditions are 
(i.) no stresses exist on the curved surface, 

(ii.) a shear Stress due to friction exists on each of the plane 
surfaces, 

and • 

(iii.) a normal compressive stress, with unknown distribution 
with respect to r, exists on each of the plane surfaces. 
Proceding from Filon's solution of displacements, the solution, 
in terms of displacements, will be assumed t o be some general function 
of r times a sine or cosine series in z plus a finite power series in 
r and z. It must also be assumed that u is odd in r and even in z and 
that w is even in r and odd in z. The assumed solutions are then 
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where f(r) must be odd in r and.g(r) must be even in r. 

Apply the displacement boundary Conditions as follows: 
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k = m v 


for m = 1,2,3*. . 


k = 
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(2.3.1) 
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Referring to equations (2.2.13) through (2.2.16), apply the 
stress boundary conditions for the cylindrical surface where r=a and 


z=z. 
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and for the plane surface where r = r and z = L/2, 
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These equations reduce to 
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Letting r go to a in equations (2.3.7) ,and (2,3.8),- 
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Eliminating p(a,^-) 
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Equation (2.3.9) is the condition at r 13 a m€-Z = L/2. If in 
equations (2;3.5) and (2.3.6), z -* L/2 the same equation results, 
therefore, the boundary cbhditions match at r = a and z => L/2. 
Consider equation (2.3.2), 
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Differentiating equation (2.3.2) with respect to S, 
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Substituting this into equation (2.3.7) 
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Consider equation (2.3.4). 
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The differential equations of equilibrium, from equations 
(2.2.11) and (2.2.12) of the previous section, became. 
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The incompressibility condition, equation (2.2.10 of the 
previous section, becomes 
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Now, there are nine unknowns, f? v (r^ , Z A^j f(r), g(r), 

^ -n — , and and there are nine algebraic equations which 
I Ar » z ' s am* 9r dz 

can be used to solve for these unknowns, equations 2.3.2, .1+, .5, .6, 

.8, .10, .12, ,13, and .14. Unfortunately, some of these equations, 
as can be easily seen, are non-linear algebraic equations and are also 
quite long. Therefore, due to the magnitude of work involved in the 
solution of these nine equations being much larger than the time 
available in which to devote to this solution, work on this aspect 
has ceased. 


2.4 Theoretical Predictions of Strains 

The strains on the surface (r = a = 1 in. ) of the cylinder, for a 
given load, will be calculated from the deformations as predicted by 
Filon's theory. The deformations are given in equations (2.1. 27) 
and (2.1.28) of section 2.1. The constants in these equations must be 
determined and can be found from equations (2.1.29) through (2.1.35). 
It has been found that Poisson's ratio for Solithane is 0.4 and in the 
following section on experimental results it has been determined from 
Table 4.1.1 that the modulus of elasticity of Solithane is 172.1 
Ib/in" ■ 

The constants have been found t o be 

A = 245.8 lb /in 2 

p = 61.46 lb/in“ 

Y = .83331 

a = 1.571 n 
n 

IjCctj) = 1.035 
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I^aj) = 1.701 

I (a) = 4.344 
1 2 

I Q (a 2 ) = 5.314 

I x (a 3 ) = 18.56 

I (a) = 20.95 
0 3 

I x ( ot^ ) = 81.46 
I (a) = 88.84 

C = 17.12, with calculations to n = 4 

f = - 11.9556 

g = - 1.573 

h = - 10,442 

Uq = .000663Q 

u : = .OOOQ16Q 

u 2 = ■ 000012Q 

D = - .OO3696Q 

E = - .0004212Q 

F - 0 00103Q 

w 0 = - .002757Q 

W 1 = .0033 m 

w 2 = - .000468Q 

D 1 = " .00187Q 

E 1 = . 000468Q 

F 1 = 0 

C = " .001069Q 
C 2 = .0008498Q 

c - - . 000139Q 
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c 4 = . 00002617Q 

A 21 = .0036T85Q 
a 22 = - .0010038Q 
A 23 = .00G20785Q 
A 24 = - .00005995Q 

A n = .005314Q 

a 12 = - .001654Q 
a 13 = .0002788Q 
A = - .00016015Q 

Using the above constants, the deformation in the radial 

direction at the center of the surface of the cylinder, u(l,0), 'Was 

found to be, with Q = T = > 

’ ^ it 

ira 

u(l»0) = 6.744 x 10“^ 

7T 

Fromthe classical theory of elasticty therefore, 

E QQ (l>0 ) = — = U(1,0). 

Using the values for applied load, Q', that were used in the experi- 
ments, the following values for shown in Table 2.4.1 can be found. 

Using the above constants, the derivative of the deformation in 
the axial direction with respect to z at r=l and z=Q was found to be > 

—4 G ' 

w'(l,0) ='- 9*802 x 10 - . From the classical theory of elasticity 

3w 

e = -r— . Using the proper values for Q r as above , e , shown in 
Z Z o z z z 

Table 2.4.2, can be found. 
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Q' (lbs) 


'00 


(x 10" ^ 
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0 
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lg. 88 

12 

25.7 6 
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38. 64 

24 

51.52 

30 

64.40 

36 

77.28 

42 

90.16 

48 

103.0k 

54 

115.92 

60 

128.80 


Circumferential Strain for Various Values of Load 

Table 2.4.1 
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91.44 

24 

121.92 

30 

152.40 
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182. 88 

42 

213. 36 

48 

243.84 

54 

274. 32 

60 

304. 80 


Axial Strain for Various Values of Load 
Table 2.4.2 
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CHAPTER 3 

EXPERIMENTAL INVESTIGATION 

3.1 Description of the Material 

Solithane resin 113 is a urethane prepolymer which upon curing 
produces solid materials. Depending upon the selection and proportion- 
ing of the catalyst employed, cured materials can be formed which vary 
from soft, rubbery compounds, through a range of intermediate states , 
to products of a hard, extremely rigid nature. 

DB castor oil was chosen as the catalyst in the preparation of 
specimens used in the experimental program here, DB castor oil is 
approximately 9 0 % triglyceride of ricinoleic acid. This catalyst wgs 
mixed with the Solithane 113 in the manner described in the following 
section and the final product was a rubbery material of intermediate 
hardness. 

Zak (10) has described the chemical aspects of Solithane resitt 
1 13 i n detail. 

3. 2 Specimen Preparation 

The repeatability of test results from one test specimen to 
another depends, to a large extent, on the modus operandi employed in 
the preparation of the specimens. Great care was taken, therefore, 
in the preparation of test specimens. The specimens utilized in this 
testing program were produced by the following method. 

The molds were of cylindrical shape, as shown in Figure 2. All 
parts of each mold were made of aluminum. The plate on the bottom of 
each mold was machined out so that the cylindrical portion of the mold 
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would fit into it. The upper plate was machined with a piston-like 
projection that would fit into the top of the cylindrical portion of 
the mold. The upper and lower end plates were attached together with 
three springs. These springs caused a force which pulled the upper 
piston- like plate down against the upper surface of the liquid. This 
was done in order to compensate for the contraction of the liquid ma- 
terial during curing. Without an end plate on the upper surface, 
shrinkage of the material during curing would have caused the surface 
of the cured material t o be somewhat concave. Also, the texture of 
the upper surface compared with the texture of the lower surface will 
differ unless a plate is used on both ends. For studies where 
friction is involved, the surface textures must be the same. 

Four molds were used for each batch of specimens prepared. The 
molds were cleaned thoroughly with a solvent. They were then sprayed 
with two mold releases, Slipicone, a product of Etw Chemical Company, 
and Fluoroglide, a Chemplast Company product. Preheating each mold 
was accomplished by placing them in an oven set at 185°F. for about 
three hours. The molds were preheated so that their temperature would 
equal that of the mixture which was going to be poured into them. 

While the molds were preheating, the basic materials used in the 
specimen preparation were mixed in the following manner. Of the total 
weight of mixture needed for one particular batch, four parts by 
weight of castor oil were vacuumed for two hours in order to remove 
as much air as possible from the liquid. The air was removed in order 
that no bubbles would form in the final solid specimens. Three parts, 
by weight, of Solith^ne 113 were added to the castor oil and the 


3 



44 


mixture of Solithane and castor oil was heated on a hot plate at a 
temperature of 130°F and simultaneously vacuumed for a period of ten 
minutes. At the end of the ten minute time lapse, a magnetically 
operated stirring device was activated and the heating, vacuuming, and 
stirring were continued for another ten minutes. At the end of the 
second ten minutes, the heat was reduced to 120°F and the vacuuming 
and stirring were continued for a third ten minute period. The object 
of the heating and vacuuming was to remove the air bubbles and the ob- 
ject of the stirring was to thoroughly mix the Solithane and castor 
oil. 

The preheated molds, without their tops, were removed from the 
oven and the mixture was poured into them to a predetermined level. 

The filled molds were then placed in a vacuum chamber and the chamber 
was evacuated for twenty minutes in order to remove air that might 
have entered the mixture during the pouring process. The tops were 
placed on the molds properly and the molds were placed into an oven 
of special design, set at 250°F. 

The design of this oven, which is shown in Figure 3, was such 
that the molds could be rotated continually during the curing time. 
This was done by attaching an electric motor, whose original 1725 rpms 
were reduced appropriately to about ten rpms, to the protrusion of a 
one-inch aluminum rod fromthe inside of the oven. The molds were af- 
fixed to this rotating rod and were therefore rotated also. This 
simultaneous rotation and heating was continued for two hours. 

The molds were removed from the oyen and their tops and bottoms 
removed. The specimens were pressed out of the cylindrical portions 
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of their molds. The mold release that had adhered to the specimens 
was removed by rubbing with a soft cloth saturated with solvent. Th$ 
specimens were then labeled appropriately as t o batch and aged for at 
least three days before being used in experiments. 

One precautionary measure was taken throughout the entire prepara- 
tion of specimens and also during storage of the final products. Since 
Solithane 113 and castor oil are very hygroscopic, the time that the 
liquids were exposed to the atmosphere, during specimen preparation, 
was reduced to a minimum. For instance, in draining the liquids from 
their storage containers, nitrogen gas was used as a forcing pressure 
in the drums. Also, the final specimens were stored at all times be- 
fore experimentation in dessicators where the dessicant used was 
calcium sulfate (Ca SO^). In addition, the experimental testing 
laboratory had a controlled atmosphere, i.e. , constant temperature at 
74°F and constant relative humidity at 45%. Therefore, any change 
that might take place in the material due to temperature or humidity 
while the experiments were being conducted was equal in magnitude. 

3.3 Experimental Apparatus and Procedure 

Because Solithane 113 is a very soft, rubber-like material, a 
method of strain measurement in which there is no physical contact 
between the specimen and the measuring device must be utilized. 

Should a device, i.e. , bonded electrical strain gages, be used where 
such physical contact occurs, the material properties of the specimen 
in the region of contact, and therefore, in the region of strain 
measurement may be altered. 
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"he strain, as yielded by the imposed testing conditions, may he. 
non-uniform, although it is assumed t o be uniform for simplicity of 
theoretical derivations of the strains on the surface. This non- 
uniformity is due to the bulging of the curved portion of the specimen 
when the axial load is applied. Local strain measurements are there- 
fore needed so that a suitable strain analysis can be performed. 

A method of measuring strains that will eliminate the two above 
mentioned experimental difficulties is the moi.r£ fringe method. The 
moire method measures strains by using the principle of mechanical 
interference of light. Two superposed grids of alternate light and 
dark lines will cause interference patterns or fringes when there is 
a relative rotation of the grids, a difference in pitch of the grids, 
a displacement between the grids or any combination of relative ro- 
tation, pitch difference, or displacement. The name, moir£, is the 
French word meaning watered- silk which describes the appearance of the 
fringes. The proper analysis, at a desired point, of the fringe 
pattern will yield the strain at that point. No physical contact is 
necessary when using the moir£ method. 

The grids used in the testing were obtained by photographic and 
diazo printing techniques. A Ronchi ruling on a one inch by two inch 
by one-quarter of an inch, piece of'glass was used as the original 
grid. Ronchi rulings are unique ija that the width of the dark lines 
are equal to the width of the light lines. The Ronchi ruling was 
used to contact print alternate light and dark lines onto Kodak, Con- 
trast Process Ortho film sheets, a high contrast, fine grain ortho- 
chromatic, antihalation film. The original grid had a pitch of 500 
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lines per inch and consequently all of the grids produced from it 
were of the same pitch. 

As noted, there are two grids employed in the moire method. One 
grid, the model grid, is affixed to the surface of the specimen t o be 
tested in the desired area of strain measurement. The second grid is 
superposed on the model grid during testing and is called the master 
grid. 

The master grid, in this case, is simply a photographic negative 
of the Ronchi ruling. A photographic negative was used as a master 
grid because its flexibility was essential when the grid was super- 
posed on the curved surface of the cylindrical specimens. 

For most materials and specimen shapes, the model grid can be 
printed directly onto the specimen surface at the desired location of 
strain analysis. Unfortunately, Solithane 113 reacts chemically with 
the solutions usually employed in printing processes. Therefore, a 
diazo printing technique had to be used, A formula for making diazo 
solution and a description of the diazo printing process can be found 
i n Appendix B. 

Since the diazo technique involves contact printing and to con- 
tact print directly onto the surface of a cylindrical specimen would 
be difficult, a bonded type gage was developed. Solithane was molded, 
by the method described previously, into a sheet 0.03 inches thick. 

The sheet was cut into smaller two inch by two inch squares. In order 
to provide contrast with the dark lines resulting fromthe diazo 
process, a ten to one mixture of white latex paint and liquid latex 
rubber was sprayed onto the surface of the small thin sheet with an 
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artist’s air brush. The liquid latex rubber was added to the latex 
paint to give the paint sufficient elasticity in order to prevent 
cracking during elongation of the model grid. Immediately after the 
paint and latex mixture was applied, the diazo solution was sprayed 
onto the sheet with the air brush. The diazo and paint- latex mixture 
was then allowed to dry, usually two to three days. Drying was done 
in a dessicator where no light could penetrate in order that the ultra- 
violet light present in the sunlight could not react with the diazo 
solution on the thin sheet. 

When dried, the coated thin sheet was placed in a vacuum printing 
frame with a photographic negative of the Ronchi ruling superposed on 
the thin sheet. A vacuum printing frame was used to ensure intimate 
contact between the thin sheet and the superposed negative during ex- 
posure. With the printing frame in operation, the diazo coating on the 
thin sheet was exposed, throught the negative, to the ultra-violet 
light source. The development process then took place as described in 
Appendix B. 

After development, the thin sheet was cut to the proper size and 
was glued to the surface of the cylindrical specimen in the area of 
desired strain measurement and with proper orientation with respect to 
the grid lineq. The glue used was General Electric #108 silicone 
cement., It is a low modulus cement as was needed because a cement 
with a modulus higher than that of the specimen would influence the 
deformation of the grid adversely. The silicone rubber cement was 
cured for at least one day before the specimen was tested. 
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Both tension and compression specimens were tested in a Tinius 
Olsen universal testing machine. This machine is shown in Figure 4. 
The maximum load available is twelve thousand pounds. Altogether, 
twelve load ranges are available for use. Tension experiments were 
conducted using the twelve pound range and compression experiments 
were run using the 120 pound range. The maximum crosshead speed is 
twenty inches per minute and is continuously variable. Both com- 
pression and tension experiments were carried out using a crosshead 
speed of one- quarter of one inch per minute. The load being applied 
to the specimen is registered on a moving strip chart that is cali- 
brated from zero to one hundred percent of the load range being 
u t i 1 ized . 

The tension specimen used was of the flat type* as shown in 
Figure 5. Its dimensions, also noted in Figure 5, conform to the 
JANAF standards. Ch the tension specimens the thin sheet with print- 
ed grid lines was glued to the test section. The thin sheet was 
glued with the grid lines perpendicular to the direction of the ap- 
plied load. The tension specimen was held in the testing machine by 
two grips like the one shown in Figure 6. The compression specimen 
used was cylindrical in shape. As shown in Figure 7, its length is 
twice its diameter. These dimensions were chosen because a specimen 
of th is size exhibits bulging when compressed but is not difficult to 
handle experimentally. 

The thin sheet was oriented on the cylindrical surface in two 
ways. For measuring axial strain the grid lines were oriented with 
the direction of the lines perpendicular to the direction of the 
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applied load. For measuring circumferential strain the grid lines 
were oriented with the direction of the lines parallel to the direction 
of the applied load. 

A cage-like apparatus, shown in Figure 4, was constructed so 
that compression could be applied to the specimen by operating the 
testing machine as though a tension experiment were being conducted. 

The following parts of the discussion on experimental procedure 
apply to tension and compression experiments alike. 

The thin sheeted area of the specimen was coated with a thin 
layer of petroleum jelly. The purpose of the jelly was to hold the 
master grid in place during the experiment. In the case of the cy- 
lindrical specimens, the master grid was held to the cylindrical sur- 
face by a girdle- like arrangement of three elastic bands. The master 
grid being affixed to the specimen, the air bubbles that had formed in 
the jelly between the grid and the thin sheet were removed by pressing 
the master grid with a soft cloth. 

The tension or compression specimen was placed in the appropriate 
grips or apparatus and the proper testing machine variables, as speci- 
fied above, were set. A Nikon F camera on a tripod, loaded with 
Kodak Tri X Pan thirty five millimeter film and equipped with exten- 
sion tubes for close up photography ,■ was positioned properly. The 
camera was used to photograph the moir€ fringes that had formed during 
experiments so that an analysis of the fringes could be performed. 

For proper illumination of the fringes during experiments, two photo- 
graphic flood lamps were used. 
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With the testing machine in operation, pictures of the fringes 
formed during deformations at specific applied loads were taken. The 
film of the fringe patterns was developed. The film negatives were 
then enlarged onto Kodak F5 or F4 Kodabromide photographic paper. 

Each enlargement was labeled properly according to test number and the 
load corresponding to the fringe pattern. 

The method of analysis of each fringe pattern is discussed in the 
following section. 

3.4 Moird Fringe Theory 

The moire fringe method of strain analysis is relatively new. Al- 
though the theory of diffraction gratings and of their manufacture date 
back to Lord Rayleigh (5) 5 1874, the application of diffraction grat- 
ings to strain measurement began with M. Dantu (2) in 1940. Since 
then, the theory involved in the analysis of the fringes formed in the 
moir4 method as well as the experimental techniques involved have ad- 
vanced rapidly. 

The formation of moire fringes is a mechanical phenomenon that 
consists of alternating light and dark regions. The fringes formed 
are not an optical phenomenon in the sense that no basic laws of 
optics are involved in their formation. Fringes result when one array 
of dots or solid lines is superposed on a second somewhat similar 
array of dots or solid lines, hereafter referred to as grids. The 
fringes are caused by a difference in the pitch of the two grids, a 
relative rotation of the grids, paralax caused by a displacement be- 
tween the grids, or any combination of a difference in pitch, a rela- 
tive rotation, or paralax. Since fringes can result from relative 
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differences in pitch and rotation, this phenomenon can be employed to 
measure strain, 

Two grids are used in the moird method of strain analysis, a 
master grid and a model grid. The master grid is a fixed array of 
straight, parallel, alternate light and dark lines with pitch, p, de- 
fined as the distance between adjacent light or dark lines of the 
grid. The model grid, a similar array of lines, is affixed to the 
specimen being tested so that as the specimen deforms under a load, 
the grid will deform with it and, therefore, will have a distorted 
pitch, p'. With the master grid superposed on the model grid, fringes 
will be formed as the model grid distorts along with the specimen. 
Although the analysis is not necessarily limited to straight, parallel 
lines, this analysis will be because straight, parallel lined grids 
were used in the experimental part of this work. 

Figure 8 shows a close up of the formation of moir6 fringes by a 
relative rotation and a difference in' pitch of two grids that are com- 
posed of straight, parallel lines. The r,s coordinate system is taken 
to be coincident with the master grid as shown. 

The acute angle between the master and model grid lines is denoted 
by 0. The angle <p is between the master grid lines and the fringe and 
is measured in the same direction as 0 and may be acute or obtuse. 

The distance between the centers of the adjacent fringes is 6. The 
pitch of the master grid is p and the pitch of the model grid is p'. 

Then, from geometrical considerations and by assuming either a 
homogeneous deformation or a sufficiently small element of a 
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nonhomogeneous deformation field, the equations for strain can be de- 
rived (4). True strain in the r-axis direction is 

OJ <3.4. l) 

*fc 

and conventional or nominal strain in the r-axis direction is 



(3.4.2) 


where 



Atf ♦ *© 


cost/) + 1 


(3.4.3) 


It should be noted that the strain determined is in the direction 
perpendicular to the direction of the grid lines of the master grid. 

As can be seen from equations (3.4.1), (3.4.2), and (3.4.3), we 
need to know 6, p, and the magnitude and direction of <j> in order to 
determine the strain and to find if the strain is compressive or ten- 
sile. In the problem in this discussion, <p = l80° and, therefore, 
equation (3.4.3) becomes 



= - 1 . 
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Therefore, equations (3.4.1) and (3.4,. 2) become 


and 





E. 


E 


(3.4.4) 

(3.4.5) 


In equations (3*4.4) and (3.4.5) > p, the pitch of the master grid 
is known and, therefore, all that must be determined in order to find 
the strain is the distance between the centers of the fringes. In 
order to find 6 from photographs of moird fringes, a graphical ap- 
proach is taken. A line is drawn, in the proper direction, across a 

£ 

series of fringes which includes the point' at which strain is to be 
measured. A plot of the accumulated fringe distances, from any con- 
venient point, versus the fringe number is made. Fringes are numbered 
consecutively with zero fringe as the point at which the accumulated 
fringe distance measurements are originated. A smooth curve is drawn 
through the plotted points. The slope of this curve at any point 
gives the value of 6, in inches per fringe, at that point, With the 
known value of 6 at any point, the value of the strain at that point 
can be found by using equations (3.4.4) and (3.4.5). 

Equations for true strain in terms of $ and the distance between 
adjacent fringe centers, measured along the coordinate axes, r and s, 
can also be derived (4), and are 


W 


eJ 


e_' 

5 

6 

r_ 


s_ 


(3.4.6) 


J 


where 
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for 


and 



1,0 < + < | 

c = (3.4.7) 

- 1 , ~ < t < TT 



(3.4.8) 



It is sometimes more accurate to measure accumulated fringe distances 
in one of the coordinate axis directions rather than to measure direct- 
ly. When this is true, equations (3.4.6), (3.4.7), and (3.4.8) are 
used to calculate strain. 

If rotation at a point is to be determined, the equations 
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if 

_ .IT 

1. 0 < <j> < 2 

c = 

, ir 

-1, - < <j> < ir 

can be used. 

Moire patterns may also be interpreted as a function of dis- 
placements as shown by Sciammarella and Durelli (7) and in an earlier 
work by J.D.C. Crisp Cl)- This interpretation of moird 1 fringes is 
based on the fact that the fringes are the loci of points presenting 
the same relative displacement in the direction normal to the master 
grid lines. 

In this method, a Cartesian system of axes, x and y, is used as 
a reference system, as shown in Figure 9* The component of displace- 
ment in the x^-axis direction is u and the component of displacement 
in the y-axis direction is V. 

A function of two variables, <f>^(x,y) } is the component of dis- 
placement of a point in a two-dimensional continuous medium parallel 
to a reference direction. The subscript, i, is 1 if the reference 
direction is the x-axis and subscript, i, is 2 if the reference 
direction is the y-axis. The function z = ij)^(x,y), in Cartesian co- 
ordinates, represents a surface. 

To determine strains from displacements, it i s necessary to com- 
pute the derivatives of the displacements. To compute the partial 
derivatives of ^(x^y), the curves of intersection of the surface 
with planes of equations X = c^ and y = c^ can be drawn, where c^ and 
c^ are constants. The procedure used to obtain a partial derivative 
is illustrated in Figure 9. The horizontal position of each point of 
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intersection of the moire' fringes with the line AB (trace of the inter- 
secting plane) is first projected on the base line CD. The orders of 
fringes are read from the moire' fringes and distances equal t o np are 
scaled up from the base line. A line drawn through the points, thus 
plotted, defines the cross-section. The slope of this curve, at a 
point, gives the derivative 
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If the preceding data are plotted directly fromthe moird pattern, 

the results obtained correspordto the Eulerian description of strain. 

For the Lagrangian description, there are two possible solutions. 

Assume that the patterns corresponding to two orthogonal directions 

have been determined. A point, p , of initial coordinates x and y 

r 0 o o 

moves, after deformation, to point p, with coordinates x and y. If 
the components of the displacements experienced by the point p^ are u 
and v, 

xo = x 1 - u 

and 

y D = y i " T • 

By using these equations, it is possible to obtain the necessary data 
to replot the moire fringes referred to the initial configuration. 

Another possibility is that the derivative of u with respect to 
can be expressed as a function of the derivative of u with respect 
to x, as follows, 

3u 3u/3lI l 

3x 1 - 3u/3x 
o 1 

• i 
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A 1 1 of the partial derivatives of the displacements with respect to the 
initial coordinates can be put into the above form. Since the partial 
derivatives of the displacements with respect to the final coordinates 
can be found, 3u/3x q , etc. can be found and the Lagrangian strains can 
therefore be found. 

A third method for analyzing strains from moir£ fringes has been 
developed recently by Sciammarella (8). This method of analysis uti- 
lizes the light intensity variations of the moir£ fringes. The new 
theory generalizes the optical law that relates the displacement field 
to the moire' fringe pattern. This yields a continuous relationship be- 
tween displacements and light intensity. Until this theory was de- 
veloped, a discontinuous relationship existed. Because of the con- 
tinuous aspect of the light intensity variation method, the precision 
of moire' fringe analysis has been increased. 

The relationship between the displacement at a point and its 
light intensity has been shown by Sciammarella t o be 


where 


u 


1 I - I ' 

arc cos 


p 


( 3 . 4 . 9 ) 


p is the master grid bitch, 

I is the light intensity at a point, 

1^ is the image average background light intensity, and 
1^ is the light intensity amplitude of the image first harmonic, 
The equipment used to measure the light intensity variations con- 
sists of a photocell mounted behind a small circular hole in the focal 
plane of a projecting system. The photocell must also be mounted on a 
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traveling carriage. A film of the fringes t o be analyzed is then in- 
serted between the photocell and a lens system. The photocell is then 
moved over the section of the fringes that are of interest and the re- 
sulting output is recorded on an X-y plotter. A sample output is shown 
in Figure 10. 

In order to compute the displacement corresponding to a point, P, 

as shown in Figure 10, it is necessary first to determine the integral 

number of fringes prior to point P. The fractional displacement is 

AB 

then determined by using equation (3.1*. 9) with I Q = I, = — and 
I = cp. The fractional displacement is then, 

Au = 27 arC COS W~~ 1 P ’ 

✓ 

The theoretical minimum displacement sensitivity for the light 
intensity variation method as calculated by Sciammarella (8) is 

At = 0.0012 p. 

The minimum strain measurable when using the geometrical approach 
or the displacement function method is a function of the gage length 
employed and the number of lines per inch on the master grid. Frojn 

i 

Figure 5 of the paper by Crisp ( 1 ^ the minimum measurable strain can 
be seen t o be approximately 0.002 inches per inch for a gage length of 
one inch and a master grid with 500 lines per inch as was used in the 
experiments in this work. 

The errors involved in the three moire strain analysis methods 
have been set down by Sciammarella. For the geometrical and displace- 
ment function methods, a bound in the error involved in strain 
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measurement when using a times 10 magnifying glass to read distances 
will be 

d£ k% . 

The error due to the measuring apparatus used in the light intensity 
variation method will be 



These errors are for measurements only and do not include errors in 
grid spacing, et cetera. 

It is apparent that the light intensity variation method is far 
superior in both sensitivity and minimum error in measurements. In 
addition, as was previously mentioned, the light intensity variation 
method yields a continuous quantity t o be measured where the two other 
methods yield discontinuous measurable quantities. However, the ex- 
pensive equipment used and the time involved in perfecting the method 
precluded the use of the light intensity variation method in this 
work. The geometrical method of moil# fringe analysis was chosen for 
use here and the light intensity variation method and the displacement 
function method were included in this discussion for informative and 
comparative reasons. 

A sample calculation of the strain from a moirl pattern can be 
found in Appendix C. 
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CHAPTER k 

EXPERIMENTAL RESULTS AND DISCUSSION 
4.1 Tension Experiments 

Tension experiments were conducted with several objectives in 
mind: (i-) to gain experience with the moire fringe method of strain 

measurement, (ii.) to verify experimental results obtained with the 
moire method with results obtained when using another method of analy- 
sis, and (iii.) to obtain a modulus of elasticity that could be used 
in the theoretical analysis. 

The third objective was reached by measuring the strain in the di- 
rection of the load application at various intervals of load. Graphs 

X, versus extension ratio, X, were plotted 
and the slope of the curve was taken as the modulus of elasticity, as 
shown in Appendix C. The moduli obtained by testing eight different 
specimens are shown in Table 4.1.1 below. As can be seen there is a 
large variation in moduli from one specimen to another, therefore, an 
average value was found. This value of these eight moduli was found 
to be 172.1 lb/in 2 and this is the value that was used in the theo- 
retical predictions of strains found, in Chapter Two, Section Four. 

The second objective was met by comparing the moduli obtained 
above with the moduli that was found by clip-gage measurements. The 
modulus obtained in the second fashion is 175.0 lb/in 2 as shown in an 
unpublished work completed here. The percent difference between the 
two is found to be 1.65%. 


of the true stress, a. 
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Specimen No. Modulus (Tb/in 2 ) 

165 E 


166 E 

167 E 

231 $ 

232 4 > 

233 $ 

236 * 

237 ♦ 


190.5 
210.0 

143.0 

162.5 

141.2 

211.0 

153.3 
165.0 


Experimentally Determined Values of Moduli 
for Different Specimens 

Table 4.1.1. 


4.2 Compression Experiments 

Uniaxial compression tests were performed on right, circular 
cylinders of moderate length which were made of cured Solithane 113 . 
The axial and circumferential strains that existed on the surfaces of 
the cylinder at the axial center were measured by means of the moir€ 
method. The range of straihs measured was determined by the charac- 
teristics of the moird grids used. The range of strains measured was 
approximately .001 in/in to .06 in/in. 

Twa types of boundary conditions were imposed on the ends of the 
cylinders. The first type imposed was such that the radial deforma- 
tions of: the end surfaces was zero. The second type imposed was such 
that the radial deformation of the end surfaces was not zero, i.e. , 
lubricant was applied to the ends in order that the friction between 
the compression plates and the ends of the cylinders was nearly zero. 
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4.3 Conclusions 

1. For relatively small compressive loads, the axial and circum- 
ferential strains are equal in magnitude for lubricated and non-lubri- 
cated ends, as shown in Figures 11 , 12, 13, and 14. 

2. With lubricated ends, the axial strain is smaller than for 
non- lubricated ends when relatively large compressive loads are applied 
to the specimens. Also, as the magnitude of the compressive load is 
increased, the difference in strains increases, as shown in Figures 13 
and 14. 

3. For a given, relatively large compressive load, the circum- 
ferential strain on a cylinder with lubricated ends is less than on a 
cylinder with non- lubricated ends, as shown in Figures 11 and 12. 

4. Strains in the axial direction are larger than strains in the 
circumferential direction, as shown in Figures 11, 12, 13, and 14. 

5. A smooth curve fromtension to compression seems to exist al- 
though further study is needed to verify this, as shown in Figure 15. 

6. The experimental results did not agree at all with the theo- 
retical results as derived by Filon (3), as shown in Figures 11 and 13. 
This was, however, expected because in the theoretical derivations the 
assumption of small strains was made and in the experimental work the 
strains were relatively large. 

7. The Poisson ratio was found to be 0.4773 for a strain rate of 
.25 in/min and at a load of 30 lbs. This agrees well with the as- 
sumption of incompressibility, where Poisson ratio is taken t o be 0.5» 
used in the derivation of the equations describing large elastic 
deformations. 
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4.4 Suggestions for Further Research 

Ch the basis of the research conducted and the results obtained 
the author makes the following suggestions for further research: 

1. That a study be made of the small strains that exist when 
small loads are applied to the specimen. 

2. That a study be made of the large strains that exist when 
large loads are applied to the specimen. 

3. That a study be made of the strains that exist at various 
points on the surface of the cylinder. 

4 . That a solution to the equations describing the strain at 
any point in the cylinder for a given load be found. 
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APPENDIX A 

MATHEMATICAL CALCULATIONS 


The equation to be solved is 


d 2 R ' . 

dR 

N 

■ 

2 . 1 

2 

— +' k 2 


2 r 

” 

2 + k 

R 2 

[dr 

dr 

l r J 

. 


Let 

R 2 = rl X + a i r + a / 2+ 


Then 



r dr 


d 2 R 

2 


2 

dr 


r y {a 0 yr -2 + a^y+ljr -1 + a 2 (y+2) + a 3 (y+3)r + 

+ a 4 (y+4)r 2 + . . . } 

r y {a Q (y-l)yr“ 2 + a i y(y+l)r -1 + a 2 (y+l)(y+2) + 
+ a 3 (y+2){y+3)r + a 4 (y+3) (y+4)r 2 + 


First let us solve the following equation. 


d 2 R 

2 

dr 


1 dR 
r dr 


k 2 R =0 


( 1 ) 


Let 

R = r y {a + a r + a r 2 + a r 3 + a r 1 * + (2) 

U 1 (L v H 

Then 

^ ^ = r y {y a Q r -2 + a ;L (^+l)r _1 + a 2 (y+2) + a 3 (y+3)r + 

+ a 4 (y+4)r 2 + . . . } . (3) 
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r U {a 0 (p-l)yr~ 2 + a 1 ii( p+l)r'" 1 + a 2 (u+l) (u+2) + a 3 (y+2) ( u+3)r 

dr 

+ a lf (U+3)(M+ i +)r 2 + (*+) 

Substitute (2), (3)> and (4) into eqn. (l) 

r lJ {a 0 ()J-l)ur- 2 + a 1 u(u+l)r" 1 + a 2 (u+l)(u+2) + a 3 ( u+2) ( u+3)r 
+ a 4 (u+3) (u+^)r 2 + . . . +u a o r ~ 2 + a ^ ( u+1 ) r" - 1 + a 2 (u+2) + 

+ a 3 (u+3)r + a tl+ (u+^)r 2 + ... " a Q r~ 2 - a^ -1 " a 2 " a 3 r 
- a 4 r 2 " ... - a Q k 2 " a^k 2 " a 2 r 2 k 2 " a^k 2 - a ft rV 2 “ 

" ...} = 0 

Equate coefficients of equal powered r 1 s , 

r " Z a Q (u-l)u + a Q u " a 0 = 0 

a Q u 2 “ a Q U + a Q u "a 0 = 0 (5) 

a Q (u 2 " 1) = 0 

r“ 1 *+ a 1 u(u+l) + a^Cu+l) - a 1 = . 0 

a i y2 + a 2 y + a i y + = 0 ( 6 ) 

a. 1 u(u+2) = 0 

r° a 2 (u+l)(u+2) + a £ (u+2) - a 2 - a Q k 2 = 0 

a 2 U 2 + a 2 3u + a 2 2 + a 2 u + a 2 2 a 2 - aj)k 2 = O 
a 2 (u+3)(u+l) - a Q k 2 =0 (7) 
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r -+ a 3 (y+2)(ja+3) + & 3 (y+3) ” a 3 " ~ 0 

a 3 y 2 + a^5y + a 3 6 + a 3 y + a 3 3 “ & 3 “ Sjk 2 = 0 
a 3 (y+4)(y+2) - ajk 2 = 0 (8) 


r 2 + a 4 (y+3)(y+4) + a^(y+4) - a 4 - a 2 k 2 = 0 

a^y 2 + a^y + a 4 12 + a 4 y + a 4 4 - a^ - E^k' 2 - 0 

a 4 (y+3)(y+5) - a 2 k 2 *= 0 (9) 

Or, looking; at (5), (6), (7)> (8), and (9) tpgether 

a Q (y 2 -l) = 0 
a 1 y ( y+2 ) = 0 
a 2 (y+3)(y+l) " a Q k 2 = 0 

a 3 (y+4)(y+2) - a^.k 2 = 0 

a 4 (y+3)(y+5) “ a 2 k 2 = 0 . 

If a 0 ^ 0, y = * 1 from equation (5); a = 0j a 3 = a^ = a 7 “ ... 

= 0, and from equation (7) 

k 2 

a 2 = a 0 (y+3) ( y+l) 

and equation (9) 

a _ a k4 

4 ° (y+3) 2 (y+l)(y+5) 


Substituting these values of a 2 > a. etc. into equation (2) 


r *■>« k?f 2 

E - r <a 0 + a 0 ( y+sTU ST 


+ a„ 


k 4 r 4 


(y+l)(y+3) 2 (y+5) 
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Now, let y = 1, a Q = k/2 

s = l kr + it* 3 * 3 + 3fe k;5r5 + 

H r = I^kr) 

The solution is invalid for y = .-l because R then blows up. 

Instead of showing how to solve for the second solution of the 
differential equation, R, = K^kr), let us prove, by using the 
Wronskian determihant, that K (kr) i s a linearly independent solution 
of the differential equation. 

It is known that, if y and y, are two solutions of a linear 
differential equation of the second order, and if yj and y'. denote 
their derivations with respect to the independent variable, then the 
solutions are linearly independent if the Wronskian determinant 

7i y 2 

y ? y' 

1 2 

does not vanish identically. 

W have that 

y , = R i = I^kr) 

y 2 = R 2 = K^kr) 
y I = k I i (kr ^ 

y ;=kK| (kr). 




TO 


Therefore, 

I L (kr) K^(kr) 

kl'(kr) kK'(kr) 

l 1 

k I (kr) K' (kr) - k IJ,(kr) K (kr) = 
k (1 1 (kr ) K^(kr)- Ij(kr) K^kr) . 


But it is known that 

I^kr) K^(kr) - I^(kr) K^(kr) = - ir- 

so that 



Then I^(kr) apd K (kr) are linearly independertt solutions of the 
original differential equations. 

Returning to the original equation to be solved, 



1 d 

+ v th* - 




( 10 ) 


£ + 1 d_ 

,2 r dr 
dr 



d 2 + 1 d_ 

_ 2 r dr 
dr 



R 


2 


= 0 


Let 


Then 


£ + 1 d_ 

, 2 r dr 
dr 





<j>(r) . 


£ + ld_ 
,2 r dr 



<J> (r ) = 0 
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where 

<j>(r) = A I (kr) + B K (kr) 

or 

4> (r ) = A I x (kr); <>(r) = B K^kr). 

Then 

1 
r 

Let 

R 2 = {a 0 + a^ + a 2 r 2 + a 3 r 3 + a^r 4 + a^r 5 + ...). 

Then 

dR 

— = air" 1 + 2a^ + 3a 0 r + la, r 2 + 5a c r 3 + ... 

r dr 2 3 4 5 

and 

d 2 R 

i. = 2a 0 + 6a r + 12a, r 2 + 20a r 3 + ... 

n o 2 3 4 5 

dr^ 

Substitute equations (12), (l3)> ( 14 ) into equation (ll) 

2a + 6a r + 12a r 2 + 20a r 3 + a r -1 + 2a„ + 3a r + la r 2 + 

23 4 51 234 

5a r 3 - a r- 2 - a r~ l - a - a r - a r 2 - a r 3 - a k 2 - a rk 2 
5 0 X ZJM-oUX 

- a 2 k 2 r 2 - a 3 k 2 r 3 - a^k 2 ** 4 - ... = A I^kr) (15) 

Ij (kr) = ^ kr + k 3 r 3 + k 5 r 5 + . . . 


(12) 

(13) 

(ll) 
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2a 2 + 6a 3 r + 12a 4 r 2 + 20a $ r 3 + a^" 1 + 2a £ + 3a $ r + 4a 4 r 2 + 5a $ r 3 
- a Q r -2 - a^ -1 - a 2 - a 3 r ~ a 4 r 2 “ a 5 r 3 - a^k 2 “ 

~ a 0 k 2 r 2 “ a 3 k' 2 r 3 -...“A^kr-A^-g k 3 r 3 " A 334 k 5 r 5 “ ... 

= 0 ‘ (16) 

Equate coefficients in equation (l6) 


r -2 

*> 

- a o 

= 0 



(17) 

r“i 

-> 

a l 

P 

I—* 

li 

O 



( 18 ) 

r° 

-> 

2a^ + 

2a 2 - 

a 2 ' 

■ a Q k 2 = 0 

(19) 



3a 2 

- a Q k 2 

= 0 





a 

2 = 0 




r 1 


6a 3 + 

3a 3 " 

a 3 

■ ajk' 2 - A i k = 0 

(20) 



8a 3 

- a^ 2 

- A 

| kr = 0 




a 3 

a i k2 

+ A 






8 



r 2 

->■ 

12a 4 

+ 4a. 
4 

“ a 

4 

- a 2 k 2 = 0 

(21) 



15a 4 

= 0 






a 4 

= 0 





r 3 20a 5 + 5a & - a g - a 3 k 2 - A ^ k 3 = 0 (22) 

24a,. = a 3 k 2 + A k 3 

t + A it k3 + A ^k 3 

a 5 ” " 24 

a k' 4 . . .3 

_ _j A k a 

a 5 192 192 


1 



73 


Place equations (17) through (22) into equation (12) 


R 2 = a^r + 


fa i k2 + Me' 


8 T l6 


r 3 + 


,a £. + Mil 

[192 192^ 


.5 + ... 


Let 


1 

2 


R 2 = - r + 


"k 2 Ak) 3 
[l 6 * iSj r + 


fk 4 A 2Ak 3 ] 

38 ? 3 §? 


r 5 + 


Let 


Then 


A = 2k 


= r + Srr kSr 3 + T§r- v 4„5 


2 2 16 


384 


k^r°+ 


or 


* - It x i (kr) ■ 


It can also be shown that 


R 2 dk K i^ kr ^ 


(23) 


(24) 


(25) 

(26) 


(27) 
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APPENDIX B 

DIAZO SOLUIION AND PRINTING PROCESS 

Initially a thin sheet was molded. Its thickness was approxi- 
mately 0.03 inches. It was made of Solithane 113 according to the 
description given in Section 3.2. The surfaces of the sheet were 
cleaned thoroughly with a solvent. 

A very thin layer of white latex paint, mixed with liquid rubber 
latex, was sprayed on one surface of the thin sheet. The rubber 
latex was mixed with the paint in order to give the paint extensibili- 
ty without cracking when the thin sheet was stretched. The purpose of 
the mixture was to provide a good contrasting background for the dark 
lines of the grid. 

Diazo printing, though not good for general photographic purposes, 
is, however, very good when reproduction of lines is needed. There- 
fore, the diazo process is readily adaptable for the application of 
grid lines to the surface of solid fuel propellants. Also, the diazo , 
process was chosen because the solutions used in more conventional 
printing processes react chemically with Solithane. 

The painted surface of the thin sheet was coated with a diazo- 
type photosensitive solution which contained the following proportion 
of chemicals (for a 100 c.c. water solution): 

Glycol . 3.0 c.c, 

Thiourea ...... 4.5 g. 

Citric Acid 4.0 g. 

Recorcinol ..... 0.3 g. 

Zinc Chloride. . . . 4.5 g. 
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Coupler III*. . . . 2.0 g. 

Diazo #8** . . . . 1.0 g. 

Diazo #9*** .... 1.0 g. 

After sufficient drying time, approximately three days, the sheet was 
placed in a vacuum printing frame with a photographic negative super- 
posed on it. The photographic negative was of the desired grid pattern. 
The printing frame was used in order to insure intimate contact be- 
tween the sheet and the negative. Illumination of the thin sheet with 
an ultra-violet light source for an appropriate amount of time destroy- 
ed the diazo salt which was not covered by the lines of the photo- 
graphic negative. Ammonia vapor was used to develop the unexposed 
portions of the diazo salt on the sheet. 


*These are trade names of the materials of the Andrews Paper and 
Chemical Co., Port Washington, N.Y. Their chemical designations 
follow. 

Coupler III is 2,3 dyhydroxy naphthalene-6*-sulfonic acid sodium 
salt. 

**Diazo #8 is para diazo NN dimethyl aniline 1/2 zinc chloride salt. 
***Diazo #9 is para diazo NN diethyl aniline zinc chloride salt. 
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APPENDIX C 

Sample Analysis of Moird Fringes 

An example of the formation of moir£ fringes is shown in Figure 
1 6 , The picture shown is a close-up of the fringes formed on a 
tension specimen, photographed when a load of 0.6 pounds had been 
applied to the specimen. The method of analysis of such a set of 
fringes is presented in section 3.4, Moire Fringe Theory. 

A plot of accumulated fringe distance yersus fringe number is 
presented in Figure 17. This plot was made, as directed in section 
3.4, from the picture of the moir£ fringes. Since the picture is a 
magnification of the original fringes, a factor of reduction must be 
included in the calculations of strains or extension ratios made from 
the data acquired from the plot of accumulated fringe distance versus 
fringe number. 

The modulus of the material tested was determined from a plot of 
true stress, a = versus extension ratio, X. In the above 

equations, P is the load applied to the specimen, A is the original 
cross-sectional area of the specimen, and h is the extension ratio. 
The slope of this plot was taken to be the modulus of the material. 
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Figure 4. Tinius Olsen Testing Machine with 
Compression Apparatus Installed. 
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Fig. 10 Sample Output of Photocell Equi] 
Light Intensity Method. 
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Fig. l4 Compressive Load Versus True Axial Strain for Lubricated Ends. 
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Figure 16. 


Moire Fringes. 
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Pig. IT Accumulated Fringe Distance Uersus Fringe Number 




